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Abstract. The status of the observations of Ultra High Energy Cosmic Rays will be reviewed, focusing on the the latest 
results of HiRes and Auger observatories. A comprehensive analytical computation scheme to compute the expected UHECR 
spectrum on earth will be presented, applying such scheme to interpret the experimental results. The phenomenological 
scenarios favored by HiRes and Auger in terms of chemical composition and spectrum will be also presented. 
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INTRODUCTION 

Ultra High Energy Cosmic Rays (UHECR) are the most 
energetic particles observed in nature with energies up 
to 3 x 10 20 eV. The detection of these particles, started 
already in the 50s with the pioneering experiments of 
Volcano Ranch in the USA and Moscow University array 
in the USSR, poses many interesting questions mainly on 
the origin and chemical composition of such fascinating 
particles. In the recent years a new step forward was done 
with the measurements performed by HiRes and AGASA 
first, and nowadays with the first five years results of the 
Pierre Auger Observatory in Argentina. 

The propagation of UHECR from the sources to the 
observer is mainly conditioned by the interaction of these 
extremely energetic particles with the intervening astro- 
physical backgrounds such as the Cosmic Microwave 
Background (CMB) and the Extraglactic Background 
Light (EBL). The study of the propagation of UHE par- 
ticles through these backgrounds is of paramount impor- 
tance to interpret the observations and to discover the as- 
trophysical origin of UHECR. Several propagation de- 
pendent features in the spectrum can be directly linked 
to the chemical composition of UHECR and/or to the 
distribution of their sources [1, 2, 3]. Among such fea- 
tures particularly important is the Greisin, Zatsepin and 
Kuzmin (GZK) suppression of the flux [1], an abrupt de- 
pletion of the observed proton spectrum due to the in- 
teraction of the UHE protons with the CMB radiation 
field. The GZK suppression, as follows from the orig- 
inal papers [1], is referred to protons and it is due to 
the photo-pion production process on the CMB radiation 
field (p + Ycmb — > 7T + p). In the case of nuclei the ex- 
pected flux also shows a suppression at the highest en- 
ergies that, depending on the nuclei specie, is due to the 
photo-disintegration process on the CMB and EBL ra- 
diation fields (A + Ycmb, ebl (A — nN) + nN). In any 
case, the interaction processes between UHE particles 



and astrophysical backgrounds will condition the end of 
the CR spectrum at the highest energies and the high en- 
ergy behavior of the flux can be used as a diagnostic 
tool for the chemical composition of the observed par- 
ticles. Another important feature in the spectrum that 
can be directly linked with the nature of the primary 
particles and on their origin (galactic/extra-galactic) is 
the pair-production dip [2]. This features is present only 
in the spectrum of UHE extragalactic protons and, as 
the GZK, is a direct consequence of the proton interac- 
tion with the CMB radiation field, in particular the dip 
brings a direct imprint of the pair production process 
p + Ycmb - ► P + e+ + e suffered by protons in their in- 
teraction with CMB radiation. 

From the experimental point of view the situation is 
far from being clear with different experiments claim- 
ing contradictory results. The HiRes experiment shows a 
proton dominated spectrum till the highest energies with 
a clear observation of the protons GZK suppression and 
the pair-production dip [4]. The chemical composition 
observed by HiRes is coherent with this picture showing 
a proton dominated spectrum at all energies E > 10 18 
eV [4]. The situation changes if the Auger results are 
taken into account. The latest release of the observed 
Auger spectrum [5] shows a not clear confirmation of 
the pair production dip and of the protons GZK suppres- 
sion. Signaling a possible deviation from a proton domi- 
nated spectrum to an heavier composition at the highest 
energies. This picture is confirmed by the chemical com- 
position observed by Auger, that shows an heavy mass 
composition at energies E > 4 x 10 18 eV [5]. 

This puzzling situation, with different experiments fa- 
voring different scenarios, shows the importance of a 
systematic study of UHECR propagation in astrophysi- 
cal backgrounds. In the present paper we will review the 
main points of the kinetic theory of UHECR propagation, 
comparing theoretical results with the latest experimen- 
tal data of Auger ad HiRes. The paper is organized as fol- 



lows: in section 2 we review the kinetic theory of propa- 
gation for protons and nuclei, in section 3 we will discuss 
different theoretical scenarios in comparison with exper- 
imental data, conclusions will take place in section 4. 

PROTONS AND NUCLEI PROPAGATION 

The propagation of charged particles (protons or nuclei) 
in astrophysical backgrounds can be suitably studied in 
the kinetic approach, that enables the computation of the 
expected fluxes on earth given any hypothesis on the 
astrophysical sources of UHECR. The main assumption 
under which the kinetic theory is build is the Continuum 
Energy Losses (CEL) approximation [6], through which 
particle interactions are treated as a continuum process 
that continuously depletes the particles energy. 

In the propagation through astrophysical backgrounds 
the interactions of particles are naturally affected by 
fluctuations, with a non-zero probability for a particle 
to travel without loosing energy. In our computation 
scheme such fluctuations are neglected; as was shown 
in [2, 7] this approach has a limited effect on the flux 
computation. Only at the highest energies (E > 100 EeV) 
fluctuations produce a deviation of the order of 10% 
respect to the flux computed with a standard Monte Carlo 
simulation [7]. 

As already anticipated in the introduction UHECR 
propagating through astrophysical backgrounds suffer 
different interaction processes loosing energy and, in the 
case of nuclei, being even destroyed. 

• protons - UHE protons interact only with the CMB 
radiation field giving rise to the two processes of 
pair production and photo-pion production. Both of 
these reactions can be treated in the CEL hypothe- 
sis. 

• nuclei - UHE nuclei interact with the CMB and 
EBL radiation fields, suffering the process of 
pair production, for which only CMB is relevant, 
and photo-disintegration, that involves both back- 
grounds CMB and EBL. While the first process can 
be treated in the CEL hypothesis, being the nuclei 
specie conserved, the second produces a change in 
the nucleus specie. Following [8] we will treat the 
photo-disintegration process as a "decaying" pro- 
cess that simply depletes the flux of the propagating 
nucleus. 

The propagation of UHE particles in the energy range 
E ~ 10 18 -i- 10 19 eV is extended over cosmological dis- 
tances with a typical path length of the order of Gpc. 
Therefore we should also take into account the adiabatic 
energy losses suffered by particles because of the cosmo- 
logical expansion of the Universe. 



The kinetic equation that describes the propagation of 
protons and nuclei are respectively: 

^E^--^[ bp (r,t)n p (r,t)} = Q p (r,t) (i) 

dn A (r,t) d n A (F,t) 

— ^[n A (r, t )b A (r,t)] + ^^ = Q A (r,t) 

(2) 

where n is the equilibrium distribution of particles, b 
are the energy losses (adiabatic expansion of the Uni- 
verse and pair/photo-pion production for protons or only 
pair-production for nuclei) Q is the injection of freshly 
accelerated particles and, in the case of nuclei, also 
the injection of secondary particles produced by photo- 
disintegration (see below). 

Concerning sources, in this paper, we will consider the 
simplest hypothesis of a uniform distribution of sources 
that inject UHECR of different species with a power law 
injection of the type: 

Gm J (r,z) = eo(z)r- r ^- r/r - (3) 

where T is the Lorentz factor of the particle, T max is the 
maximum Lorentz factor that the sources can provide 
and Yg is the power law index of the distribution of the 
accelerated particles at the sources. In our discussion, de- 
pending on the particle specie, we will consider different 
injection parameters Qq and T max , keeping the power law 
index y g fixed for all species. 

The energy losses b for protons or nuclei depend only 
on the CMB field and in the CEL hypothesis can be 
computed analytically through [2]: 

b{Y) = ±f de r a(e r )f(e r )e r f de?™^l (4) 
2r j £o Je r /2r e z 

where T is the Lorentz factor of the UHE particle, e r is 
the energy of the background photon in the rest frame of 
the UHE particle, £o is the threshold of the considered 
reaction in the rest system of the UHE particle, <j(s r ) is 
the cross-section of the process, f(e r ) is the mean frac- 
tion of energy lost by UHE particle in a single collision 
in the laboratory system, and n CMB {e) is the density of 
the CMB background photons in the Laboratory refer- 
ence frame. 

The energy losses (4) do not take into account the ex- 
pansion of the Universe and are intended at the present 
epoch (z = 0), in order to generalize energy losses at any 
red-shift z, because in the computation of b only CMB 
plays a role, one can simply use the following recipe 
b(r,z) = (1 +z) 2 b((l +z)r,z = 0). The expansion of 
the universe is taken into account by adding b(F 7 z) + 

r(l +z)H(z), with H(z) =/W(l +z) 3 £} m +£iA. In the 
present paper we will always assume standard cosmol- 
ogy with H = 72 Km/(s Mpc), Q. m = 0.27 and Q. A = 
0.73. 
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FIGURE 1. protons path length, with labels of the relevant 
process that fix the behavior in different energy regimes. The 
sharp suppression of the path length responsible of the expected 
GZK cut-off is clearly seen as a consequence of the photo-pion 
production process. Only the CMB radiation field is relevant 
for proton propagation (see the text). 



In figure 1 we show the behavior of the path-length 
A = c/b„ (at z = 0) for protons, from the figure it is 
evident the effect of the photo-pion production process 
with a sharp decrease of the contributing Universe that 
is at the basis of the expected GZK cut-off in the proton 
spectrum. 

In the case of nuclei the energy losses term bA in 
equation (2) depends only on the pair-production process 
on CMB and on the universe expansion. The rate b^ can 
be computed as in equation (4) and it is easy to show that 
it is related to the protons pair-production loss rate by the 
relation: 

b A (T) = \b b {Y) (5) 

where Z,A are the charge (atomic number) and atomic 
mass number of the nucleus. The second process that 
affects nuclei propagation is photo-disintegration over 
CMB and EBL backgrounds. This process is treated as a 
decaying process that depletes the flux of nuclei, it enters 
in the kinetic equation (see equation (2)) through a sort of 
"life-time" of the nucleus under the photo-disintegration 
process. This "life-time" corresponds to the mean time 
needed to a nucleus of Lorentz factor F and atomic 
number A to lose, at least, one of its nucleons: 

— = -^j[ de r a{e r ,A)v{e r )e r [ ^ £ 'W £ ) 

T A 2T 2 J £o (A) Je r /(2r) £ 2 

(6) 

where a{e ri A) is the photo-disintegration cross-section 
and v(e r ) is the molteplicity associated to this process, 
namely the average number of nucleons extracted from 
the nucleus by a single interaction and tibkg = ucmb + 
kebl- The dependence on red-shift of %a directly fol- 
lows from the evolution with red-shift of the background 
photon densities iicmb and uebl- In the case of CMB 
this dependence is known analytically while for the EBL 




FIGURE 2. Nuclei path length under the processes of pair 
production X = c/b^ (blue dotted) and photo-disintegration 
A = c%A' In the case of photo-disintegration we have considered 
the path-length computed taking into account all relevant back- 
grounds CMB and EBL (red continuos line) and only CMB 
(black dotted line). The size of the Universe is also shown as 
greed dashed line. 



we will refer to the models presented by Stecker [9]. 
In figure 2 we show the behavior of the path length A 
for different nuclei species, being A = c/bA, for pair- 
production, or A = cTa, for photo-disintegration. In fig- 
ure 2 we show also the effect of the EBL by plotting A 
computed taking into account both backgrounds (contin- 
uos line) and only CMB (black dotted line). It is impor- 
tant to note that the photo-disintegration process starts 
to contribute to the propagation of nuclei at a Lorentz 
factor that is almost independent of the nuclei specie 
r cr ~ 2 x 10 9 [8]. This is an important general character- 
istic of nuclei photo-disintegration process from which 
we can immediately deduce the dependence on the nu- 
cleus specie of the energy corresponding to the photo- 
disintegration suppression of the nuclei flux: 



Ef m = Am N T cr 



(7) 



being A the atomic number of the nucleus and the 
proton mass. From equation (7) it is evident how the 
flux behavior could bring informations on the chemical 
composition of the UHECR, in the case of Helium (A = 
4) the suppression is expected around energies E ~ 10 19 
eV while in the case of Iron (A = 56) the suppression is 
expected at higher energies E ~ 10 20 eV. 

Let us discuss now the generation function g^(r,f) 
in the rhs of Eq. (2). One should distinguish among pri- 
mary nuclei, i.e. nuclei accelerated at the source and in- 
jected in the intergalactic space, and secondary nuclei 
and protons, i.e. particles produced as secondaries in 
the photo-disintegration chain. In the case of primaries 
the injection function has the form (3), while the injec- 
tion of secondaries should be modeled taking into ac- 
count the characteristics of the photo-disintegration pro- 
cess. The dominant process of photo-disintegration is the 



one nucleon (AO emission, namely the process (A + 1 ) + 
Jbkg A+N, this follows directly from the behavior of 
the photo-disintegration cross-section (see [8] and ref- 
erences therein) that shows the giant dipole resonance 
corresponding to one nucleon emission. Moreover, at the 
typical energies ofUHECR(E> 10 17 eV) one can safely 
neglect the nucleus recoil so that photo-disintegration 
will conserve the Lorentz factor of the particles. There- 
fore the production rate of secondary A— nucleus and 
A— associated nucleons will be given by 



G.(r,z) = ^(r, z ) = ^£) 

' z A +i[T,z) 



(8) 



where T^+i is the photo-disintegration life-time of the 
nucleus father (A + 1) and n A +\ is its equilibrium distri- 
bution, solution of the kinetic equation (2). 

Using equation (8) we can build a system of cou- 
pled differential equations that starting from primary in- 
jected nuclei (Ao), with an injection given by (3), follows 
the complete photo-disintegration chain for all secondary 
nuclei (A < Ao) and nucleons. Clearly secondary pro- 
tons 1 propagation will be described by the proper kinetic 
equation (1) with an injection term given by (8). The so- 
lution of the kinetic equation for protons and nuclei can 
be worked out analytically. In the case of protons: 



n p (T,i 



(1 



(9) 



being Q p the injection of primary protons (equation (3)) 
or secondary protons (equation (8)) and F' = r'(r,z) is 
the characteristic function of the kinetic equation [8]. In 
the case of nuclei: 



«A(r,z) 



Q A (T' z )— e -^ r '^ 

-z')H{z'r Ay ,Z) dT 



(10) 

being, again, Q A the injection of primary nuclei (3) or 
secondary (8). The exponential term in Eq. (10) is given 
by 



e -*u(rV) =exp 



dz" 



1 



{\+z")H{z") T A (r",z") 



(11) 

in which one easily recognizes the survival probability 
during the propagation time t ' — t for a nucleus with fixed 
A. Therefore, exp(— 77) provides the suppression of large 
z' in the integral in Eq. (10). The derivative term dY' /dY 
present in both solutions (9) and (10) can be computed 
according to [8]. Finally, the maximum red-shift of inte- 
gration Zmax can be directly computed solving the equa- 
tion r'(r,z) = E max , assuming a maximum acceleration 




1 Neutrons decay very fast into protons, so we will always refer to 
secondary protons. 



FIGURE 3. Pure proton spectrum in comparison with HiRes 
data [4], data are from HiRes II monocular (boxes) and HiRes 
I monocular (triangles). The two curves r/ fof and f] ee show 
the total spectrum and the spectrum calculated including only 
adiabatic and pair-production energy losses. 



energy E max — AmNY max (A = 1 for protons) in the injec- 
tion of freshly accelerated particles. 



COMPARISON WITH OBSERVATIONS 

Using the solutions of the kinetic equation (9) and (10) 
we can compute the expected flux on earth (at z = 0). 
In this section we will discuss the two alternative sce- 
narios that follows from the observations of Auger and 
HiRes. Let us first consider the case of HiRes data. In 
figure 3 we show the comparison of these data with a 
pure proton spectrum, the comparison is performed by 
means of the modification factor 77(E) [2]. This quantity 
is given by the ratio of the energy spectrum n p (E) calcu- 
lated from equation (9) (with all energy losses taken into 
account), and the unmodified spectrum n unm (E), where 
only adiabatic energy loss (Universe expansion) is in- 
cluded: 77(E) = n p {E)/n unm {E) [2]. The modification 
factor is a convenient quantity to describe, in a model in- 
dependent way [2], any feature in the spectrum related to 
energy losses. From figure 3 it is evident the agreement 
of the HiRes observed spectrum with the expected pair 
production dip, typical of a pure proton spectrum with a 
best fit injection power law y g — 2.7 [2]. 

The nature of the spectrum suppression seen in figure 
3 can be further characterized taking the integral spec- 
trum as observed by HiRes. The GZK suppression for 
protons in the integral spectrum can be tagged by the en- 
ergy scale E\n, that corresponds to the energy at which 
the integral spectrum J(> E) becomes half of the sim- 



pie power-law extrapolation spectrum KE~Y. Theoretical 
calculations for protons predict a value E\/2 = 10 19,72 eV 
for a wide range of injection power laws y g = 2.1 -j- 2.8 
[10]. The measured value of E\ji as claimed by the 
HiRes collaboration is in fairly good agreement with 
the predicted value, being E^,f es = l0 19 - 73±007 . From 
these observations of HiRes, with some caution, one may 
conclude that HiRes has detected the two main signa- 
tures of a pure proton composition: the pair production 
dip and the GZK suppression. This important conclusion 
of the HiRes observations is confirmed from the direct 
measurement of the UHECR chemical composition per- 
formed by this experiment. In particular, the HiRes mea- 
surement of the elongation rate X max and its root-mean- 
square (RMS) confirm a strong dominance of protons in 
the UHECR composition at all energies E > 10 18 eV [4]. 

The picture emerging from HiRes observations is not 
confirmed by Auger. The Auger data on mass composi- 
tion and spectra are quite different from those of HiRes. 
In particular, the Auger observations show a nuclei dom- 
inated spectrum already at energies around 4 x 10 18 eV. 
This conclusion is quite robust in the Auger data and it 
follows from both observations on elongation rate as well 
as on its RMS. The change with energy of mass composi- 
tion in the Auger data is quite smooth with a progressive 
increase of the heavy nuclei content starting already at 
2 x 10 18 eV. 

Coherently with the elongation rate observations the 
UHECR flux observed by Auger is difficult to recon- 
cile with a pure proton spectrum. In figure 4 we show 
the comparison of the Auger spectrum with a pure pro- 
ton spectrum obtained fixing the injection power law 
Y g =2.8 and with different choices of the maximum at- 
tainable energy at the source (as labeled). It is clear that 
the observed spectrum is well described by a pure proton 
flux only at energies below (3-^5) x 10 18 eV. At larger 
energies the theoretical proton spectrum is not compati- 
ble with the Auger observations [11], that don't show the 
typical dip behavior as discussed in the case of HiRes 
data. The same conclusion holds changing the injection 
power law, taking a very flat injection with y g = 2.0 the 
situation doesn't change with the theoretical proton spec- 
trum compatible with Auger observations only at low en- 
ergy, E< (3 -T- 5) x 10 18 [1 1]. 

As discussed in the case of HiRes observations, an- 
other important signature of the proton content at the 
highest energies is related to the shape of the flux sup- 
pression. Also under this respect the flux observed by 
Auger seems not compatible with a pure proton compo- 
sition. The value of E^n and the behavior of the sup- 
pression are better related to an heavy nuclei composition 
than a pure proton. 

As recently pointed out in [11], the experimental ev- 
idence emerging from Auger data can be understood in 
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FIGURE 4. Pure proton spectrum with an injection power 
law Yg = 2.8 and for different choices of the maximum energy 
Emax (as labelled) in comparison with the observed Auger flux. 



terms of different primary species and their maximum 
energy at the source. A large set of astrophysical acceler- 
ation mechanisms, in particular those based on diffusive 
shocks, show a maximum energy at the source which is 
rigidity dependent, being related to the electric charge of 
the accelerated particle: T^ lax °< ZF max with Tf nax max- 
imum Lorentz factor of a nucleus with atomic number 
Z and Tmax maximum Lorentz factor for protons. Us- 
ing this characteristic of the maximum energy it is pos- 
sible to explain the progressively heavier composition 
with increasing energy observed by Auger. This goal is 
reached by the so-called "disappointing model" [11], that 
describes the Auger observations by means of a multi- 
component injection of UHECR at the sources. At en- 
ergy higher than ZEf nax nuclei with charge Z' < Z dis- 
appear, while heavier nuclei with larger Z survive. Start- 
ing from Emax — (4-^10) xlO 18 eV, the higher energies 
are accessible only for nuclei with progressively larger 
values of Z. In particular, the maximum observed en- 
ergy must correspond to Iron nuclei, which can reach 
E? n e ax ~ (1 3) x 10 20 eV. In figure 5 we plot the UHECR 
spectrum in a two component model: protons and Iron, 
in comparison with the Auger observed flux, for an in- 
jection power law index y g = 2.0 and a value of the max- 
imum energy for protons of Emax = 4 x 10 18 eV. From 
this figure it is evident the good agreement of this simple 
(two-components) model with the Auger observations. 
The predictions of this model are very disappointing for 
future detectors. In fact, the maximum energy for Iron 
Emax — (1 x 10 20 eV implies an energy per nucleon 
below the threshold E N < E F m % x jA ~ (2 5) x 10 18 eV, 
i.e. well below the GZK energy. Therefore, the produc- 
tion of secondary detectable particles such as neutrinos 
or gamma-rays will be subdominant respect to the photo- 
disintegration process of nuclei and no alternative detec- 
tion of the UHECR propagation will be possible [11]. 
Moreover, it will be impossible any correlation study of 
the UHECR events with possible astrophysical sources 
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FIGURE 5. UHECR spectrum in a two component model 
with protons and Iron nuclei, the injection power law index is 
y g = 2.0 and the maximum energy is E max = 4Z x 10 18 eV. 



eV. This scenario emerging from Auger observations, if 
confirmed, will be quite disappointing because the dom- 
inance of nuclei at the highest energies will seriously 
harm any experimental study of correlation with sources 
as well as any detection of UHE neutrinos or gamma-rays 
produced by UHECR propagation. 

Let us conclude by stating that the experimental ob- 
servation of the UHECR chemical composition at the 
highest energies has a paramount importance in choos- 
ing among the two alternative scenarios depicted in this 
paper and establishing the future directions of this field 
of research. Unfortunately the available observations at 
the highest energies are still affected by poor statistics 
and renewed experimental efforts are needed in order to 
unveil the nature of UHECR. 



because the effect of galactic magnetic field will sub- 
stantially deviate the nuclei trajectories at the highest en- 
ergies [11]. For these reasons the model at hand, while 
gives a good description of the Auger observations, was 
named "disappointing" by the authors of [1 1]. 

CONCLUSIONS 

In conclusion, we face at present the most serious dis- 
agreement in the observational data of the two biggest 
experiments in UHECR. HiRes observes signatures of 
proton propagation through CMB in the form of the pair- 
production dip and GZK cutoff. Moreover, these obser- 
vations are well confirmed by the HiRes direct measure- 
ment of a proton dominated mass composition. The study 
of UHECR propagation enables us to accommodate 
HiRes observations in a pure proton model, with sources 
characterized by a steep injection spectrum y„ = 2.5-^2.8 
and an high maximum energy Em ax = (1-^3) x 10 20 eV. 

On the other hand, Auger clearly observes an high- 
energy steepening of the spectrum, but its position and 
shape are rather different from the prediction of the GZK 
cutoff. Moreover, the behavior of the spectrum observed 
by Auger in the energy range 1 x 10 18 -f4x 10 19 eV 
doesn't confirm the pair-production dip typical of pro- 
tons, signaling a substantial nuclei contamination in the 
flux observations. Coherently with this, the mass com- 
position directly observed by Auger at E > 4 x 10 18 eV 
shows a dominance of nuclei that becomes progressively 
heavier increasing the energy and reaches a pure Iron 
composition at E ~ 10 19 eV. 

From our analysis follows that the Auger data favor a 
multi-component spectrum at the sources with a conser- 
vative explanation in terms of flat injection y g = 2.0 -f- 2.3 
and a relatively low maximum energy for protons Ef nax ~ 
(3-j-5) x 10 18 eV, that corresponds to a maximum en- 
ergy for Iron nuclei of the order of E^f ax ~ (1 -=-2) X 10 20 
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